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The one-dimensional Bose gas is an unusual superfluid. In contrast to higher spatial dimensions, the existence 
of non-classical rotational inertia is not directly linked to the dissipationless motion of infinitesimal impurities. 
Recently, experimental tests with ultracold atoms have begun and quantitative predictions for the drag force 
experienced by moving obstacles have become available. This topical review discusses the drag force obtained 
from linear response theory in relation to Landau's criterion of superfluidity. Based upon improved analytical 
and numerical understanding of the dynamical structure factor, results for different obstacle potentials are ob- 
tained, including single impurities, optical lattices and random potentials generated from speckle patterns. The 
dynamical breakdown of superfluidity in random potentials is discussed in relation to Anderson localization and 
the predicted superfluid-insulator transition in these systems. 
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I. INTRODUCTION 

Superfluidity in a neutral gas or liquid is not easily defined 
but rather understood as a complex cluster of phenomena. The 
associated properties may include frictionless flow through 
thin capillaries, suppression of the classical inertial moment, 
metastable cuiTents, quantized circulation (vortices), Joseph- 
son effect (coherent tunneling), and so on (see, e.g., Refs. [1- 
6]). There is a close and deep analogy between superfluidity 
in a neutral system and superconductivity in a charged system 
[4, 6]. 

The 3D weakly-interacting Bose gas has all the super- 
fluid properties mentioned above, which can be inferred from 
the existence of an order parameter represented by the wave 
function of the Bose-Einstein condensate (BEC). By contrast, 
there is no BEC in a repulsive ID Bose gas even at zero tem- 
perature in the thermodynamic limit, provided interactions are 
independent of particle velocities [7, 8]. This can be easily 
proved using the Bogoliubov theorem [7]. This pre- 

dicts a l/q^ divergence at small momentum in the average 
occupation number rig for nonzero temperature and 1/q di- 
vergence for zero temperature. Nevertheless, the existence of 
BEC is neither a sufficient nor necessary condition for super- 
fluidity [1,6], and a one-dimensional system of bosons may be 
superfluid under some conditions. However, whether a system 
is superfluid or not depends very much on how superfluidity 
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is defined, because one-dimensional systems may exhibit only 
some but not all of the superfluid phenomena. 

Here we study superfluidity in an atomic gas of repulsive 
spinless bosons in the ID regime of very narrow ring confine- 
ment. The investigations focus mainly on the metastability 
of the circulating-current states in various regimes. However, 
we also discuss another important aspect of superfluidity rel- 
evant to a ID system, which is the non-classical moment of 
inertia or Hess-Fairbank effect [9] and the quantization of cir- 
culation. As we argue below, a perfect Hess-Fairbank effect 
and quantization of circulation occur for the homogeneous gas 
of repulsive spinless bosons in one dimension, while metasta- 
bihty of currents does not, in general. Note that the Hess- 
Fairbank effect is much easier to investigate than metastability 
of current because of its "equilibrium" nature [6, 10]. Indeed, 
it can be explained with the properties of the low-lying en- 
ergy excitation spectrum of the system due to the abiUty of 
the system to relax to the ground state in the reference frame 
where the walls (i.e. the trapping potentials) are at rest (see 
Sec. II). The metastability of currents is a much more com- 
plicated phenomenon to study, because at sufficientiy large 
gas velocities, the system is obviously not in the ground state 
but in a metastable state. In order to study such an effect, one 
needs to understand transitions between states, which presents 
a more intricate problem. 

Ideally, in order to study the decay of ring currents in a con- 
trolled manner, the gas should be kept in a ring or a torus-like 
geometry with defined defects. The defects may cause tran- 
sitions to the states of lower energies, thus leading to energy 
dissipation, related to a friction force; this is called the drag 
force. The question of metastability then becomes equivalent 
to the drag force of a small and heavy impurity that is dragged 
through the resting gas. In the scope of this topical review, 
we consider mainly infinitesimal impurities and calculate the 
lowest order terms in linear response of the interacting gas to 
perturbation by the impurity. The authors of [1 1] followed a 
different approach by calculating the effects of finite impuri- 
ties on the flow of a weakly-interacting Bose-Einstein conden- 
sate. 

In spite of rapid progress in experimental techniques along 
this line [12-14], so far no conclusive experimental data on 
the drag force or metastability of ring currents in the ID Bose 
gas is available, and thus, this is one of the outstanding funda- 
mental questions remaining about the properties of ultra-cold 
Bose gases [2]. Not long ago, an experiment along this Une 
was carried out [15], in which the propagation of spin im- 
purity atoms through a strongly interacting one-dimensional 
Bose gas was observed in a cigar-shaped geometry. The mo- 
tion of the center-of-mass position of the wave packet is de- 
scribed fairly well by the drag force, calculated with the dy- 
namic structure factor of the Bose gas in the regime of infinite 
boson interactions. In the recent experiment of Ref. [16], the 
dynamics of light impurities in a bath of bosonic atoms was 
investigated and the decay of breathing mode oscillations was 
observed. In another fine of experiments, atoms were sub- 
jected to a moving optical lattice potential and the momen- 
tum transfer was measured [17-19]. This implies that one 
can experimentally obtain the drag force of a specific exter- 



nal potential acting on the gas. Experiments were also done 
with ultra-cold atoms in random and pseudo-random poten- 
tials. The direct observation of Anderson localization was re- 
ported in Refs. [20-22]. In particular, spreading of a ID Bose 
gas in artificially created random potentials was experimen- 
tally investigated [20]. Below we show that the superfluid- 
insulator phase diagram of such a system can be obtained by 
calculating the drag force. 

The notion of drag force turns out to be theoretically fun- 
damental, because it generahzes Landau's famous criterion of 
superfluidity. According to Landau, an obstacle in a gas, mov- 
ing with velocity v, may cause transitions from the ground 
state of the gas to excited states lying on the fine e = pv in 
the energy-momentum space. If all the spectrum is above this 
line, the motion cannot excite the system, and it is thus su- 
perfluid. However, it is also possible that even when the line 
intersects the spectrum, the transition probabilities to these 
states are strongly suppressed due to boson interactions or to 
the specific kind of external perturbing potential. In this case, 
the drag force gives us a quantitative measure of superfluidity. 

This paper is organized as follows. The basic model of the 
ID Bose gas considered in this paper is introduced in Sec. II. 
In the subsequent section, we study the Hess-Fairbank effect 
and its relation to the Landau criterion of superfluidity. In 
Sec. Ill, we derive an expression for the drag force through 
the dynamic structure factor within linear response theory and 
show that the notion of the drag force generalizes the Lan- 
dau criterion of superfluidity. The values of the drag force 
in various Umiting regimes are calculated in Sec. IV. In the 
subsequent section, the Luttinger hquid theory is exploited to 
describe the drag force at small impurity velocities. An ex- 
act method for obtaining the drag force for a finite number of 
bosons, exploiting the algebraic Bethe ansatz, is briefly dis- 
cussed in Sec. V C. In Sec. V D we consider a simple interpo- 
lation formula for the DSF, which works for arbitrary strength 
of the interparticle interactions. In the subsequent section, 
we show that another approach to the one-dimensional Bose 
gas, the theory of phase slip transitions [23], reproduces the 
same power-law behaviour for the drag force as the Luttinger 
liquid theory does. In Sec. VF, the drag force in the limit 
of infinite interactions is obtained analytically beyond Hnear 
response theory. The decay of ring currents is examined in 
Sec. VI. In Sec. VII A, the zero-temperature phase diagram 
is obtained for the superfluid-insulator transition of the one- 
dimensional Bose gas in moving shallow lattices for arbitrary 
values of velocity, filling factor, and strength of boson inter- 
actions. The motion of the ID Bose gas in random potentials 
and the related phase diagram at zero temperature are consid- 
ered in Sec. VII B. Finally, in the conclusion the results and 
prospects are briefly discussed. 



n. LANDAU CRITEMON OF SUPERFLUIDITY AND 
HESS-FAIRBANK EFFECT 

Model - We model cold bosonic atoms in a waveguide-like 

micro trap by a simple ID gas of bosons with point inter- 
actions of strength qb > (i-e. the Lieb-Liniger (LL) model 
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[24]) 

^ ^2 Q2 

and impose periodic boundary conditions on the wave func- 
tions. The strength of interactions can be measured in terms 
of the dimensionless parameter 7 = mg^/ {h?n), where n is 
the linear density and m is the mass. In the limit of large 7, 
the model is known as the Tonks-Girardeau (TG) gas [25, 26]. 
In this limit, it can be mapped onto an ideal Fermi gas since 
infinite contact repulsions emulate the Pauli principle. In the 
opposite limit of small 7, we recover the Bogoliubov model 
of weakly interacting bosons. For an overview of theoretical 
approaches to the one-dimensional Bose gas see the recent re- 
view [27]. 

Landau criterion - In the LL model the total momentum is a 
good quantum number [28], and periodic boundary conditions 
quantize it in units of 2Trh/L, where L is the ring circumfer- 
ence. Classification of all the excitations can be done [29] 
in the same manner as for the Tonks-Girardeau gas of non- 
interacting fermions (7 +00). Thus, in order to create an 
elementary particle-like excitation, one needs to add a quasi- 
momentum beyond the occupied Fermi segment. By contrast, 
for a hole-like excitation, one needs to remove a quasimomen- 
tum lying inside the Fermi segment. All the excitations can 
be constructed from the above elementary excitations. Due to 
the conservation of total number of particles, the number of 
particle-like excitations coincides with that of hole-like exci- 
tations. The low-lying spectrum of iV = nL bosons as shown 
in Fig. 1 has local minima [30] at the supercurrent states / 
(/ — 0,1,2,...) with momenta pi = 2TrnhI and excitation 
energies 

Si=pj/{2Nm). (2) 

These correspond to Galilean transformations of the ground 
state with velocities vj = pj/{Nm). The minima do not 
depend on interactions and tend to zero in the limit of large 
system size at constant density. The first supercurrent state 
is also called the umklapp excitation [24] by analogy with 
periodic lattices because it can be reached from the ground 
state by imparting the momentum HKy to each particle, where 
= 27r/_L is the reciprocal wave vector in the ring geom- 
etry. As explained above, this changes the total momentum 
while preserving the internal state of the system. 

Suppose that the gas is put into rotation with linear veloc- 
ity vi, and after that, is braked with an artificial macroscopic 
"obstacle," e.g., created by a laser beam [31]. In the refer- 
ence frame where the gas is at rest, the obstacle moves with 
velocity vj. In a superfluid we expect to see no energy dis- 
sipation, and the drag force is zero (the current is persistent). 
Otherwise one can observe decay of the current. It follows 
from energy conservation that the transitions from the ground 
state caused by the moving obstacle with velocity 1;, lie on 
the line uj — vk in the energy-momentum plane. According 
to Landau, if the excitation spectrum lies above this line, the 
motion cannot excite the system, which is then superfluid. As 



is seen from Fig. 1, the Landau critical velocity (when the 
line touches the spectrum) equals Vc — vi/2 = h'K/{mL). 
This implies that any supercurrent state with J ^ 1 is unstable 
since vi > Vc- However, in 3D we also have similar supercur- 
rent states, which apparently leads to the paradoxical absence 
of current metastability. The solution to this is that we need 
to consider not only the spectrum but also probabilities of ex- 
citations. Below we argue that in the 3D case, the probability 
to excite supercurrents is vanishingly small, while in the ID 
case it depends on the strength of bosonic interactions. 




FIG. 1. (Color online) Schematic of the excitation spectrum of the 
ID Bose gas in a perfectly isotropic ring. The supercurrent states 
/ lie on the parabola h^k'^ /{2Nm) (dotted line). Excitations oc- 
cur in the shaded area; the discrete structure of the spectrum is not 
shown for simplicity. The blue (dark) area represents single particle- 
hole excitations [29]. Motion of the impurity with respect to the gas 
causes transitions from the ground state to the states lying on the 
straight (red) line. Reproduced from Ref. [32], Copyright (c)2009 
The American Physical Society 

Hess-Fairbank effect - When the walls of a toroidal con- 
tainer are set into rotation adiabatically with a small tangen- 
tial velocity vd, a superfluid stays at rest while a normal fluid 
follows the container. This effect leads to a nonclassical ro- 
tational inertia of superfluid systems, which can be used to 
determine the superfluid fraction [6, 10]. For the ID Bose 
gas, rotation of the annular trap amounts to shifting the exci- 
tation spectrum to e — vdp as shown in Fig. 2. It is assumed 
that an unspecified relaxation mechanism allows the system 
to relax to the ground state in the frame where the trap is at 
rest. The low-lying LL excitation spectrum is a convex func- 
tion of momentum for ^ p ^ pi [29], and, hence, the mo- 
mentum zero state remains the ground state for |wd| < Vc. 
This leads to the Hess-Fairbank effect for the ID Bose gas 
for arbitrary repulsive interactions 7 > [33, 34]. According 
to this equilibrium property which is completely determined 
by the low-lying energy spectrum [1], the ID Bose gas has 
a 100% superfluid fraction and zero rotational inertia at zero 
temperature. The same results were obtained by using instan- 
ton techniques [35] and within Luttinger liquid theory for a 
finite number of bosons [34, 36]. 

It is the convexity of the low-lying excitation spectrum 
between the supercurrent states (see Fig. 2) that allows us 
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to obtain this result without numerical calculations. The 
minima of energy can be reached only in the supercurrent 
states (see Fig. 2b), whose energy are known analytically 
with Eq. (2). In the usual way (see, e.g., Ref. [34]), the 
superfluid component can be defined and obtained numer- 
ically through the second derivative of the phase in the 
twisted boundary conditions ■ ■ ■ ,Xj + L, . . . , x^) = 

exp{i(p)ip{xi, . . . ,Xj, . . . ,xn), j = 1 A^, imposed on 
the wave functions. We emphasize that this procedure is 
equivalent to the method used in this paper, because the 
Galilean transformation implies a phase gradient of the wave 
functions observed in the frame of the moving walls, and the 
gradient leads to twisted boundary conditions. It is also con- 
sistent with the definition of the superfluid component in three 
dimensions through the Fourier-transformed autocorrelation 
function of the transverse current [6, 10]. 

Note that the Bose-Fermi mapping for the TG gas implies 
that periodic boundary conditions for the Bose gas only trans- 
late into periodic boundary conditions for the Fermi gas for 
odd N but into antiperiodic ones for even TV [37], in con- 
trast to the true Fermi gas with periodic conditions. One can 
see that for even N, there is no Hess-Fairbank effect in the 
true Fermi gas (although the TG gas always shows the Hess- 
Fairbank effect). This is due to instability of fermions at the 
Fermi point, which results from the degeneracy of the ground 
state for an even number of fermions. 



III. DRAG FORCE AS A GENERALIZATION OF 
LANDAU'S CRITERION OF SUPERFLUIDITY 

In order to study frictionless motion, let us consider an im- 
purity of mass TOj, moving with velocity vina medium of par- 
ticles. For the major part of this topical review, we will treat 
the interaction between the impurity and Bose gas perturba- 
tively in linear response theory and thus reduce the problem to 
calculating properties of the integrable ID Bose gas, although 
in general the system will lose integrability when adding an 
impurity. Using Fermi's golden rule, one can easily show (see 
App. A) that the resulting friction leads to an energy loss per 
unit time given by 
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(27r 
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N 



(3) 



Here Vi{q) is the Fourier transform of the interaction potential 
Vi{r) between the impurity and the particles, D is the spa- 
tial dimension, n is the density of particles (number of parti- 
cles per D-dimensional volume), and S{q, uj) is the dynamical 
structure factor (DSF) of the medium. It is given by the defi- 
nition [5] 

S{q,uj) = ^e-^^"|(m|,5pq|n)|25(fia; - + E„,), 

(4) 

with Z = J2m 6xp(— /?i?m) being the partition function and 
/3 being the inverse temperature. Here 5pq = e"'^ ''^ — 
A (q) is the Fourier component of the operator of the density 
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FIG. 2. (Color online) (a) Low-lying excitation spectrum emin(fc) for 
= 10 particles for ID repulsive bosons. At 7 <C 1, the low-lying 
spectrum verges towards that of the ideal Bose gas, which lies on the 
straight segments between points / = 1, 2, 3, . . .. (b) Quantization 
of current velocity for ID repulsive bosons under influence of a mov- 
ing trap. Shown are the low-energy excitations of the ID Bose gas in 
the moving frame £min(A:) — vohk, calculated from the Bethe-ansatz 
equations [24] for different values of the coupling strength (compare 
with emin(fc) of Fig. 1). Inset: The velocity of the gas at equilibrium 
changes abruptly at values of driving velocity vo/vc = 1, 3, 5, . . ., 
since the gas occupies the state with lowest energy. In particular, the 
system is at rest when the driving velocity is less than = v-p /N 
(perfect Hess-Fairbank effect). Here, fep = nn, vy = hkY/m, 
and vf = h'^kp/{2m). Reproduced from Ref. [38], Copyright © 
Siberian Federal University. 



fluctuations, A(q) = 1 at q = and A(q) = otherwise. At 
zero temperature, the structure factor takes a simpler form 

^(<7,w) = ^|(0|(5/5q|n)|25(ftc^-S„ + So). (5) 



The DSF relates to the time-dependent density coiTelator 
through the Fourier transformation 



S{q,u)^N 



dtd^r 



{Spir,t)Sp{0,0))/n, 



(6) 
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where Sp{r, t) = S{r — rj{t)) — n is the operator of the 
density fluctuations. The DSF obeys the /-sum rule [5] 
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duiijS{q,uj) = Nq^/{2m). 



(7) 



The drag force is defined by the formula E = —F^ ■ v. In 
this paper, we will use the expression (3) for a heavy impurity 
V ^ hq/m.[ in one dimension. It yields for the drag force 
_Fv = dqq\Vi{q)\'^S{q,qv)/L. Using the properties of 

the DSF S{q, w) = 5'(-g, w) = el^'"^S{q, -w), which follow 
from its definition (4), we obtain 
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dqq \mq)\^S{q,qv)[l-eM-0J^<lv)]/L. (8) 



This is the most general form of the drag force within lin- 
ear response theory. The form of the impurity interaction po- 
tential can be of importance as, e.g., in the cases of shallow 
lattices and random potentials discussed in Sees. VllA and 
VII B, respectively. When the impurity interaction is of short- 
range type, we can replace it with good accuracy by a contact 
interaction Vi{r) = giS{x), which leads to Vi{q) = Qi and 
yields at zero temperature [39] 



_2gfnm _ , 



dqqS{q,qv)/L. (9) 



In the first equaUty, we introduce the dimensionless drag force 
/v with the help of the "natural" unit 2gfnrn/ ti^. lis physical 
nature will be discussed below in Sees. IV A and V F. 

The notion of drag force generalizes the Landau criterion 
of superfluidity. Indeed, the integral in Eq. (9) is taken along 
the line uj = qv in the uj-q plane. If the excitation spectrum 
lies completely above the line then the integral vanishes, as 
one can see from the DSF definition (5). On the other hand, 
the integral can be infinitesimaUy small or vanish even if the 
spectrum lies below the line but the excitation probabilities, 
given by the corresponding matrix elements Vi{q){0\6pq\n), 
are suppressed. The drag force can vanish even if the system 
is not superfluid in principle, but excitations are not accessed 
with a given potential. This happens when the matrix ele- 
ments {0\6pq\n) are not small but the Fourier transform of the 

impurity potential Vi{q) takes non-zero values only in a finite 
region of q-space. We consider such an interesting case in 
Sees. VII A and VllB below. 

Thus, determining the drag force within linear response the- 
ory is reduced to the problem of calculating the dynamic struc- 
ture factor. Below we consider various approximations for the 
dynamic structure factor and the associated drag force. 



rV. THE DRAG FORCE IN DIFFERENT REGIMES 

A. Large impurity velocities 

Let us consider the case of contact interactions with the 
impurity, Vi{r) = giS{x). At large impurity velocity v ^ 



fmn/ra, the main contribution to the integral in Eq. (9) comes 
from the high momentum region of the DSF, which can be 
calculated analyticaUy [40]. Indeed, at large velocities, the 
momentum transfer from the impurity to the particles is big 
enough to neglect the interparticle interactions. Then one can 
use the DSF values of the ideal gas [5, 40] 



np6 ( huj 



h^q-^ h^pq 
"2m 
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(10) 



Here Up = {a'pttp) is the average occupation numbers of par- 
ticles, and the plus is taken for bosons and the minus for 
fermions. The second equaUty in Eq. (10) is due to the large 
momentum hq ~ mv » hnn, which leads to UpUp+q ~ for 
arbitrary values of p. Substituting Eq. (10) into Eq. (9) yields 
the value of the drag force 2g?mn/h^, which can be used as 
the force unit. 

It is not difficult to see that this result is vaUd beyond Un- 
ear response theory, given by Eq. (9). Indeed, at sufficiently 
large particle velocity, the initial particle momentum can be 
neglected during the scattering. Then in the reference frame 
where the impurity remains at rest, the relative particle mo- 
mentum is hq ~ mv. Therefore, the reflection coefficient is 
determined by the squared absolute value of the scattering am- 
plitude in the Bom approximation m^gf /h^q"^. Each particle, 
being reflected, transfers momentum 2hq. The total number 
of scattered particles per unit time is nv = nhq/m. The prod- 
uct of the last three quantities yields the value of momentum 
transfer per unit time, that is, the drag force. This quantity, in- 
dependent of the wave vector, is nothing else but the limiting 
value of the drag force, obtained above. 



B. The Tonks-Girardeau regime 

Let us first investigate the drag force for the TG gas (7 — >■ 
+00), having the same structure factor as the ideal Fermi gas. 
The DSF at zero temperature is well-known in the thermody- 
namic limit [41, 42] 



S{k,cj) 



£F 

N 



(11) 



fora;_(fc) < w < oj+{k), and zero otherwise. Herea;±(A;) = 
h\2kYk ± fc^|/(2m) are the limiting dispersions that bound 
quasiparticle-quasihole excitations (see Fig. 4). By definition, 
= 7rn and = h^k^/{2m) are the Fermi wave vector and 
energy of the TG gas, respectively. As follows from Eq. (11) 
for the DSF, the transition probability from the ground state is 
inversely proportional to the momentum transfer but does not 
depend on the excitation energies within the borders a;±(fc). 
Simple integration in Eq. (9) then yields 



/v = 



u/i;f, ^ w ^ Up, 



1, 



(12) 
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where vp = hirn/m is the sound velocity in the TG regime. 
The result is represented in Fig. 3a. The TG gas is obviously 
not superfluid. 




FIG. 3. (Color online) The dimensionless drag force (9) versus im- 
purity velocity in various approximations, c is the speed of sound, 
(a) Tonks-Girardeau and Bogoliubov regimes, (b) The random phase 
approximations (RPA) at non-zero temperatures. The curve becomes 
smoother when the temperature grows. 



C. The drag force in the Bogoliubov regime 

The opposite regime to the TG gas is the limit of weak in- 
teractions 7^1 (the Bogoliubov regime). The crucial point 
in Bogoliubov theory [43] is the developed Bose-Einstein con- 
densate. In spite of the absence of Bose-Einstein condensation 
in one dimension [8, 44], the upper dispersion curve w+(fc) is 
well described at small 7 [29] by the Bogoliubov relation 



huji. 



r2 + 4^fceF7/T^ 



(13) 



where Tk = h'^k'^ /{2m) is the free-particle energy spectrum. 
The Bogoliubov theory also yields the correct values of the 
ground state energy and chemical potential. This paradox can 
be explained [45] by a strong singularity of the DSF near a;+ 
in the Bogoliubov regime. As a result, it is localized almost 
completely within a small vicinity of the upper branch. This 
tendency can be seen even at 7 = 10 (Fig. 4). Thus, the 
behavior of the DSF simulates the (5-function spike, which ap- 
pears due to the Bose-Einstein condensate. One can simply 



put S'Bog(fc, w) = C6{u! — LUk) and determine the constant C 
from the /-sum rule (7) 



r{k,uj) = N- — S{uj - LOk). 



(14) 



The drag force is then obtained analytically: it is a step 
function 



/v = 



0, sC w s; c, 

1, w ^ c, 



(15) 



shown in Fig. 3a. The sound velocity in the Bogoliubov 
regime is given by c = wf^/t/tt. At non-zero temperatures 
the result does not change. This result was first obtained in 
Ref. [39] by means of the GP equation. Note that the super- 
sonic impurity motion in three dimensions was earlier studied 
by the same method in Ref. [46]. 



D. The linear approximation near the Tonks-Girardeau 
regime 

For finite 7, the model can also be mapped onto a Fermi 
gas [37] with local interactions, inversely proportional to 
[41, 42, 47, 48]. Using the explicit form of the interactions, 
one can develop the time-dependent Hartree-Fock scheme 
[41, 42] in the strong-coupling regime with small parameter 
1/7. Approximations of the linear response functions on this 
level are known as Random Phase approximation (RPA) with 
exchange or generalized RPA [4, 49]. The scheme yields the 
correct expansion of the DSF up to the first order [41, 42] 



N 4:k \ 7 / 27 w| 



(16) 

for a;_(fc) ^ w ^ uj+{k), and zero elsewhere. The symbol 
0{x) denotes terms of order x or smaller The limiting dis- 
persions in the strong-coupling regime take the form 



2m 



V 7/ \7 
From the linear part of the dispersion (17) at small momentum 
uj±{k) ~ ck, we obtain the well-known result [29, 34] for 
the sound velocity at zero temperature c = fF(l — 4/7) + 
0(1/7'). 

In the same manner as in Sec. IV B, we derive from Eqs. (9) 
and (16) 

,ln(i;/c)" 



1 



7 



1, 



c 

V ^ c. 



(18) 



As expected, the linear approximation (18) works badly for 
small values of impurity velocity. This is due to the anoma- 
lous behaviour of the DSF within the linear approximation in 
vicinity of the umklapp excitation point at w = 0, fc = 2fcF. 
Indeed, Eq. (18) leads to unphysical negative values of the 
drag force at sufficiently small impurity velocities. For this 
reason, we need a more careful examination of the drag force 
in this regime. 
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V. THEORETICAL APPROACHES 

A. Random phase approximation near tlie Tonlcs-Girardeau 
regime 

The DSF in the RPA was calculated and described in details 
in Sec. IVB of Ref. [42]. The RPA is based on the Hartree- 
Fock scheme, which is appropriate at sufficiently large inter- 
actions between bosons 7 > 8. The equation for the drag 
force within the RPA is derived in App. B. 

The RPA approximation has some advantages. It always 
gives positive values of the drag force and is applicable at 
non-zero temperatures. The results are shown in Fig. 3b. The 
RPA scheme can be extended to non-homogeneous systems 
beyond the local density approximation [50]. However, the 
RPA scheme also works badly in the vicinity of the umk- 
lapp excitation. As a consequence, the formula for the drag 
force (B7), obtained within the RPA, does not reproduce the 
correct power-law behaviour at small velocities, see discus- 
sions in Sees. V B and V D below. Another disadvantage of 
the scheme is that the drag force as a function of velocity has 
an unphysical peak near f = c at zero temperatures, which 
becomes quite pronounced for 7 < 15. This is an apparent 
artifact of RPA. In the next sections we consider much bet- 
ter approximations that work for all values of the interaction 
strength and impurity velocities. 



B. Luttinger liquid tlieory 

Luttinger liquid theory allows us to correctly describe the 
low-energy excitations of a ID system of particles whose in- 
teractions are independent of the velocities. It faithfully han- 
dles the nonlinear effects of these interactions and yields val- 
ues of the DSF at low energies for arbitrary coupling strength 
[30, 39] and thus can be used to calculate the drag force at 
small values of the impurity velocity. 

The behaviour of the DSF in the vicinity of the umklapp 
excitation is related to the long-range asymptotics of the time- 
dependent density-density correlator. The asymptotics for 
k-px ^ 1 and t ^ x/c 



{dpix,t)dp{0,0)) 



K 



+ ■ 



1 



Air'^rfl \{x- ctY {x + ctY 

cos{2kYx) 



A,{K)- 



'{x^-cH^)K 



(19) 



can be justified by generahzing Haldane's method [30] or us- 
ing perturbation theory [51, 52] or quantum inverse scattering 
method [53] or conformal field theory [54, 55]. The first two 
terms are related to the behaviour of the DSF in vicinity of 
w = 0, g = 0, while the third term is related to the umk- 
lapp excitations. Substituting the third term of Eq. (19) into 
Eq. (6) yields the DSF in the vicinity of the "umklapp" point 
(fc = 2kF = 27vn, w = 0) [39] 



5'(fc,a;) 
N 



nc 



huj 



2K 



M{K) 1 



ujl{k) 



K-l 



for u! > uj^{k), and zero otherwise. Within Luttinger- 
liquid theory, the dispersion is linear near the umklapp point: 
(jj-{k) ~ c|fc — 27rn|. By definition, K = h-Kn/{mc) is the 
Luttinger parameter. For repulsive bosons, the value of pa- 
rameter K lies between 1 (TG gas) and -|-oo (ideal Bose gas). 
The value of the Luttinger parameter in the strong-coupling 
regime 



if = 1 + 4/7 + 0(1/72) 



(21) 



is derived with the expression for the sound velocity obtained 
in Sec. IV D. The coefficients in Eqs. (19) and (20) are related 
by [56] 



2K 



MK) 



(22) 



with T{K) being the gamma function. The coefficient Ai {K) 
is non-universal, and is explicitly known for the Lieb-Liniger 
model [57, 58]. We wiU make use of its value in two lim- 
iting cases: Ai{K) = st K = 1 and Ai{K) ~ 
4i-3^exp(-27cis:)7r/r2(if) for if » 1 [39, 56], where 
7c = 0.5772 ... is the Euler constant. 

By comparing the first-order expansion (16) in the vicinity 
of the umklapp point with Eq. (20) and using the expansion 
(21), one can easily obtain [45] the model-dependent coeffi- 
cient at large but finite interactions when if — 1 <C 1 



Ai{K) 



1- (1-F41n2)(if - 1) +0{{K-lf). 



(23) 



Using the DSF (20) predicted by Luttinger theory and equa- 
tion (9) for the drag force, we arrive at the expression for the 
drag force at small velocities <C c [39] 



T{K) 



T{K + 1/2) 



(20) 



In the TG regime, this formula yields the same values of the 
drag force as Eq. (12). Equation (24) gives us the universal 
exponent of the power-law behaviour of the drag force at small 
velocities: ^ v"^^^^. The same result was obtained in 
Ref. [36]. While the non-universal coefficient ^i(if ) is now 
known for arbitrary strength of interactions [57-59], its actual 
expression is too unwieldy to be considered here. In Sec. V D 
we prefer to consider another approach, which allows us to 
determine it approximately. 



C. Tlie algebraic Bethe ansatz and ABACUS 

The exact integrabihty of the Lieb-Liniger model now per- 
mits the direct numerical calculation of dynamical correlation 
functions such as the DSF [60] for systems with finite num- 
bers of particles using the ABACUS algorithm [61]. The strat- 
egy consists in using the Lehmann representation (5). The 
eigenstates themselves are constructed by the Bethe Ansatz. 
The state norms [62, 63] and matrix element of the density op- 
erator [64, 65] being exactly known, the sum over intermedi- 
ate states can be taken starting from the dominant few-particle 
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0.5 1 1.5 2 2.5 3 

k/kp 



FIG. 4. (Color online) Numerical values of the DSF (4) for the 
coupling parameter 7 = 10 [60]. The dimensionless value of the 
rescaled DSF S{k,uj)eF/N is shown in shades of gray between zero 
(white) and 1.0 (black). The upper and lower solid (blue) lines repre- 
sent the dispersions uj+{k) and lu- (fc), respectively, limiting the sin- 
gle "particle-hole" excitations in the Lieb-Liniger model at T = 0. 
The dispersions are obtained numerically by solving the system of 
integral equations [29]. The gray scale plot of the DSF demonstrates 
that the main contribution to the DSF comes from the single particle- 
hole excitations, lying inside the region aj_ (fc) ^ lj ^ uj+ (k) (see 
also Fig. 1). Only one point at fc = kc, shown in full (red) cir- 
cle, contributes to the integral when the perturber is a shallow cosine 
potential with a reciprocal vector kc- Reproduced from Ref. [32], 
Copyright (5)2009 The American Physical Society 



CUS data for the DSF. The derivative with respect to v can 
then be taken numerically by fitting an interpolating polyno- 
mial to the data points for i^^ (t;) in the vicinity of the velocity 
for which the drag force needs to be calculated. The resulting 
data for the drag force are illustrated in Fig. 5a. 

The advantage of this method is that the reliability of the re- 
sults is more or less independent of the value of the interaction 
parameter, in the sense that states and matrix elements can be 
individually constructed irrespective of what 7 is. The remain- 
ing issues are the distribution of coiTelation weight among the 
excitations, and how this affects speed of convergence. The 
limit of very small 7 is the easiest to treat, since only very few 
states in the vicinity of the Type 1 mode are of importance. 
The DSF then becomes almost a (5— function in the Type 1 
dispersion relation, and the drag force tends to the step func- 
tion as expected in this limit. The opposite case of infinite 7 is 
also straightforward, since then only single particle-hole ex- 
citations have non-negligible matrix elements, the drag force 
becomes a constant between Type 1 and 2 dispersions, and the 
drag force acquires a constant slope. For interaction values be- 
tween these two extremes, many eigenstates must be summed 
over for good convergence, and the DSF takes on a nontriv- 
ial lineshape making the drag force take the positive-curvature 
shapes in Fig. 5a. 



D. An effective approximation for the dynamic structure 
factor and drag force 



states and let run until a satisfactory saturation of sum rules is 
achieved. The results for the DSF take the form illustrated in 
Fig. 4. For generic values of 7, most of the signal is con- 
centrated between the dispersions uj+ (fc) and a;_ (fc) of single 
particle/hole (Type 1/2 in Lieb's notation) excitations. The 
signal is by construction identically zero below w_ (fc) (since 
no eigenstates are found there); above w+(fc), the signal is 
caiTied by multiple particle-hole excitations and is very weak 
in view of the small matrix elements of the density operator 
between these states and the ground state. 

We use data for iV = 150 particles (7 = 5 and 20), N = 
200 (7 = 1) and N = 300 (7 = 0.25). The /-sum rule 
saturations at fc = 2fcF were 99.64% (7 = 20), 97.81% (7 = 
5), 99.06% (7 = 1) and 99.08% (7 = 0.25), the saturation 
getting better at smaller momenta. The drag force is computed 
from the numerical DSF data in the following way. Since the 
DSF in finite size is given by discrete but densely distributed 
i5-function peaks, we consider the integral 

2 poo pqv 



HL 



i- dqJ2mSp,W)\^QifiQV-En + Eo), (25) 

JO ^ 



whose derivative with respect to v simply gives the drag force 
(9). Here 9 is the Heaviside step function. The integral in 
ui conveniently gets rid of all energy (5-functions in the ex- 
pression for the DSF, and this integrated quantity is readily 
computed without need for smoothing using the raw ABA- 



In Ref. [45, 66], an interpolating expression was suggested 
for the DSF; 



^(fc,^) = C 



(w" -UJ°^)^'- 



(26) 



fora;_(fc) ^ w ^ a;+(fc), and S{k,uj) — otherwise. Here, 
l^+{k) and fJ.-{k) are the exact exponents [67][68] 



S{k,L 



- w±(fc) 



(27) 



at the borders of the spectrum w+(fc) and w_(fc). We also 
put by definition a = 1 + The most general way 

of obtaining w±(fc), ^±{k), and K is to solve numerically 
the corresponding integral equations of Refs. [29] and [67], 
respectively. 

It follows from energy and momentum conservation that 
S{k,Lu) is exactly equal to zero below uj^{k) for ^ fc ^ 
27m. In the other regions of w > a;+ and u < uj^ (for 
fc > 2TTii), possible contributions can arise due to coupling 
to multi-particle excitations [29]. However, these contribu- 
tions are known to vanish in the Tonks-Girardeau (7 — > 00) 
and Bogoliubov (7 — > 0) limits and are found to be very small 
numerically for finite interactions [60]. 

The exponents /x± are non-negative [69]. As a conse- 
quence, the DSF diverges at the upper branch cj^. At the 
lower branch a;_, the DSF shows a continuous transition to 
zero for any finite value of 7 except for the specific point 
7 = +00 (or K = 1) of the Tonks-Girardeau gas, where the 
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FIG. 5. (Color online) The dimensionless drag force versus the ve- 
locity (relative to the sound velocity) of the impurity at various values 
of the coupling parameter 7. (a) The solid (blue) lines represent the 
force obtained with Eqs. (9) and (26), open circles are the numerical 
data obtained using ABACUS [60]. (b) The values of the dimen- 
sional drag force obtained from the interpolation formula (26) [solid 
(blue) lines] are compared with that of the analytical formula (29) for 
small velocities [dashed (red) lines]. 



DSF remains finite but has a discontinuous transition to zero 
at both boundaries a;_ and lo^. 

The normalization constant C depends on momentum but 
not on frequency and is determined from the /-sum rule (7). 
The expression (26) is applicable for all ranges of the param- 
eters fc, Lj, and 7 with increasing accuracy at large 7 [45]. 

The parameter a is needed to reconcile the limiting value 
of the exponent ^_ (27m) — 2\/K{^/K — 1) in the vicinity of 
the umklapp point [67] and the Luttinger theory predictions, 
given by Eq. (20). Now one can see from (26) that 



S{k,uj) 



^2(if-l)^ 



k = 2?™, 
k ^ 2Tm. 



(28) 



Thus, the suggested formula (26) is consistent with both the 
Luttinger liquid behavior at the umklapp point and Imam- 
bekov's and Glazman's power-law behavior in its vicinity, as 
it should be. A more detailed discussion can be found in 
Ref. [45]. Similar approximations for the DSF of ID Bose 
gas, confined in a harmonic trap, are considered in Ref. [70]. 

The drag force can now be calculated by means of Eqs. (9) 
and (26) for arbitrary strength of interactions and arbitrary ve- 
locities. The results are shown in Fig. 5a. 

For the important question whether persistent currents may 
exist at all, the small velocity regime is most relevant, which 



is dominated by transitions near the first supercurrent state 
(umklapp point at lu — and k = 2kp). The drag force 
in this regime has a power-law dependence on the velocity 
Fv ~ w^^^^ for ti ^ c, as first found by Astrakharchik and 
Pitaevskii [39]. From Eqs. (9) and (26) we obtain 



/v 



2K 



(i) 



2K-1 



4er 



huj-i-(2kp) 



2K 



r(i + ^ - M+(2fcF)) r(i + /i_(2fcF))r(i + ^) 
Y(^)r(i-M+(2fcF)) r(i + M_(2fcF) + ^) 



(29) 

where T{x) is Euler's Gamma-function, and /i_(2fcF) = 
2\/K{\/K — 1) [67]. This formula is valid for arbitrary cou- 
pling constant and works even in the Bogoliubov regime at 
7 ^ 1. In practice, Eq. (29) works well up to u < 0.1c, see 
Fig. 5b. The strong suppression of the drag force in the Bo- 
goliubov regime appears because of the large exponent value 
2K— 1 3> 1 and the large argument 2K/ a ^ 1 of the gamma 
function in the denominator. Certainly, drag force values of 
order 10^^*^ lie beyond realizable experiments and as a matter 
of fact, tell us about superfluidity in this regime. 



E. Drag force from the phase slip transitions 

The quantum number /, describing the supercuiTent with 
energy (2), is nothing else but the phase winding number of 



the original Bose field 



'ne 



[30]. The transitions be- 



tween the supercurrent states are the phase slip transitions, 
changing the phase winding number. One can write down ex- 
plicitly the low-energy Hamiltonian [23] that consists of the 
supercurrent states and phonon excitations, interacting due to 
a static impurity. The interaction Vi{x) between the impurity 
and bosons can be rewritten in terms of the supercurrent states 
and the phonons. One can then directly derive the transition 
probability per unit time between the supercurrent states / and 
/' from Fermi's golden rule [71]. The result reads at zero tem- 
perature [see Eq. (8) of Ref. [71]] 



ftr[l + a]r[i + a] iir 



^ , (30) 



where a = {I — I') K and ^/// = £/ — £//. The transition 
probability (30) is determined up to an overall factor, because 
the dimensionless parameters g and 7, which are of the order 
of 1, cannot be defined exactly from the long-range effective 
Hamiltonian. The energy eq is a characteristic high-energy 
cutoff for the phonon spectrum. 

For the /th supercuiTent decay at zero temperature due to 
the static impurity, the dominant contribution comes from the 
transition to the (/ — l)th state, and one can write for the 
energy loss per unit time 



\99in 



hT\l 



K]T[\ + K] V 7^0 



;/,/-! 



2K 



(31) 
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For large winding number /, we have the relation ^7.7- 1 ~ 
2TTnhvj, resulted from Eq. (2). By putting the energy cutoff 
eo to be equal to 2'Knhc and remembering the definition of the 
drag force \E\ — F^v, we obtain from Eq. (31) 



F„ 



\ggin\ 



4^3/2 



1 



2K 



(32) 



Thus, we come to the dependence Fy ^ v"^^ ^ at small ve- 
locities, in agreement with the previous results (24) and (29). 



F. Direct calculation for the Tonks Girardeau gas beyond 
linear response theory 

In this section we consider the TG gas, or the ideal gas of 
spin-polarized fermions. In order to find the drag force, we 
choose a frame of reference, in which the massive impurity 
is at rest. Moving particles are then scattered by the impurity 
potential Vi{x) — gi6{x). For an ideal Fermi gas the scatter- 
ing process for each particle is independent of the scattering 
of the other particles. The reflection coefficient for a particle 
with wave vector k is easily determined from the one-body 
Schrodinger equation: r{k) = m^gf /{m^gf + h*k'^). If dA''^ 
is the number of particles with wave vectors lying between 
k and k + dfc, then these particles transfer momentum to the 
impurity per unit time dF = 2h^nr{k)k'^dNk/m. Since in 
the chosen frame of reference, the momentum distribution is 
the Fermi distribution shifted by mv, we can write down the 
explicit expression of the drag force 



2r^ 

m 



dfc 
2^ 



n{k — mv/h) sign(fc)r(fc)fc^ 



(33) 



where n{k) = 1 for — vrn ^ fc ^ 7m, and zero otherwise. The 
integral is readily taken 



2mng? 



7r(1^4) 



^ + ^ (^arctan 
1 ^ 



arctan 



^ (arctan ^'■^ — arctan ^'^'^^^^ 



(34) 



Here we denote 71 = gim/ (fi^n), and ^ = v/vp. The results 
are shown in Fig. 6. As discussed in Sec. IV A, the limiting 
value of the drag force at large velocities coincides with that 
obtained within linear response theory. At small impurity cou- 
pling 7i <C 1, Eq. (34) reproduces the linear response formula 
(12), see Fig. 3a. 



VI. CONSQUENCE OF DRAG FORCE: VELOCITY 
DAMPING 



In the presence of an obstacle as, e.g., shown in Fig. 7, a 
ring current can decay due to successive transitions to super- 
current states with smaller momentum. Starting in one of the 




FIG. 6. The diagram shows the dimensionless drag force (9) for 
the TG gas at various values of the impurity coupling parameter 
7i = gim/{K^n). The results are obtained beyond the linear re- 
sponse approximation. Note that the absolute value of the drag force 
Fv is proportional to 7^ . 




FIG. 7. An experimental scheme to observe current decays. An ar- 
tificial "impurity" is created by optical methods. It can be switched 
on adiabatically or abruptly. 



local minima of the excitation spectrum as seen in Fig. 1, the 
kinetic energy of the center-of-mass translation will be trans- 
formed into elementary excitations above a lower supercur- 
rent state while still conserving the total energy. The elemen- 
tary excitations are quasiparticle-quasihole excitations in the 
Bethe-ansatz wave function [29]. Both in energy and char- 
acter most of these excitation lie between the phonon-like 
branch and the w_ branch, which is related to dark soli- 
tons (see Fig. 4). Assuming that these excitations have lit- 
tle effect on successive transitions, we estimate [32] the de- 
cay of the center-of-mass velocity v by the classical equa- 
tion Nmv = ~Fy{v), where F^ is given by Eqs. (9) and 
(26). This equation was integrated numerically and the result 
is shown in Fig. 8. At the initial supersonic velocity, where 
the drag force is saturated (see Fig. 5a) the supercurrent ex- 
periences constant deceleration. For v < c the drag force 
decreases and consequently the deceleration slows down. For 
the TG gas we find an analytical solution for exponential de- 
cay v{t) = woexp(— i/r) for vq ^ vp. In the weakly- 
interacting regime, the decay may be slow compared to ex- 
perimental time scales. 
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FIG. 8. (Color online) Decay of the ring current velocity of ID 
bosons from the initial velocity of l.lvp at t = 0. The solid (blue) 
and dashed (red) lines represent the results obtained with the ap- 
proximate formula and ABACUS, respectively. The time scale is 
r = NTvh^/{2mgf). Reproduced from Ref. [32], Copyright (c)2009 
The American Physical Society 



Vn. DRAG FORCE IN EXTENDED POTENTIALS 
A. ID bosons in a moving shallow lattice 

Equation (8) can be verified experimentally for different 
types of obstacles: for Vi{x) — giS{x) all the points at the 
transition line contribute to the drag force, while for the pe- 
riodic potential with the spatial period a only a set of dis- 
crete points in the uj-k plane do [32, 38]. This is simply due 
to the fact that a periodic potential has only a discrete set of 
Fourier components with momenta jkQ with j being integer 
and fc(3 = 27r/a being the reciprocal lattice vector. For in- 
stance, we have two nonzero components kc and — fcg for 
optical lattice potential Vi{x) — gi, cos{2TTx/a) 



mq)\'/L - ngliSiq - kc) + d{q + kc)) /2 



(35) 



in the thermodynamic limit (n = const, L — 7> oo). The fill- 
ing factor of the lattice, that is, the number of particles per 
site, is given by a = 2TTn/kG, because the total number of 
lattice periods equals Lkc/i^Tr). Substituting Eq. (35) into 
the general expression for the drag force (8), we obtain at zero 
temperature 



F^^7:glkGS{kG,kGv)/2. 



(36) 



The values of the obtained drag force can easily be estimated 
with the ui-k diagram for the DSF, see Fig. 4. Note that the 
drag force is now proportional to the total number of parti- 
cles, by contrast to the point-like impurity case, described by 
Eq. (9). This is because the lattice potential is non-local. 

Let us recall that the model considered in this paper as- 
sumes periodic boundary conditions. However, according to 
the general principles of statistical mechanics, the boundary 
conditions do not play a role in the thermodynamic limit. 
Hence, the obtained formula (36) can be exploited at suffi- 
ciently large number of particles even in the case of a cigar- 
shaped quasi-lD gas of bosons. This equation gives us the 




FIG. 9. (Color online) Zero temperature phase diagram for 
superfluid-isolator transition of the Bose gas in a moving shallow 
lattice: dimensionless drag force F^fi? / {gl^mL) versus the lattice 
velocity (in units up) and the interaction strength 7 at various values 
of the filling factor a. The dimensionless values are represented in 
shades of gray between zero (white) and 1.0 (black). The solid (blue) 
lines correspond to the DSF borders u}+(k) and uj- (k), respectively. 



momentum transfer per unit time from a moving shallow lat- 
tice, which can be measured experimentally [17-19]. 

In the case of Bogoliubov and TG regimes, the drag force 
admits analytic solutions. As discussed in Sec. IV C, at small 
7 the upper dispersion curve Lo+{k) is described well by the 
Bogoliubov relation (13), and the behavior of the DSF sim- 
ulates the (5-function spike in accordance with Eq. (14). We 
derive from Eqs. (14) and (36) 



g^mL 1 
h? 2ai 



1 



7 



(37) 



where by definition ^ = v/v-p- 

In the TG regime, discussed in Sec. IV B, the DSF is given 
by Eq. (11) with the limiting dispersions ljj± being known an- 
alytically. Then Eqs. (11) and (36) yield 



(38) 
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where we put by definition ^-t = 1 1 ± l/a| . 

The values of the drag force (36) obtained from the inter- 
polating expression for DSF (26) is shown in Figs. 9 and 10. 
The filling factor a — 1 (kc — 2kp) corresponds to the Mott 
insulator state in a deep lattice. As can be seen from Fig. 4, 
at this value of the reciprocal vector the DSF is almost inde- 
pendent of u when 7^1, and the drag force takes non-zero 
values for arbitrary v ^ aj+ {kG)/kG, see Fig. 9. 

By contrast, at small 7 its non-zero values practically local- 
ize in the vicinity of t; = Lo+{kG)/kG, as shown in Fig. 10. 
Then superfluidity breaks down when the point cj = vkG and 
k = kG lies exactly on the Bogoliubov dispersion curve, see 
Fig. 4. Taking into account that the Bogoliubov dispersion 
is very close to the free particle one, we obtain for the break 
point fcc?^' = hkQ/{2m). One can see that this point coincides 
with the point of dynamical instability for bosons in the cosine 
shallow lattice. Indeed, the dynamical instability appears at a 
value of lattice velocity v corresponding to a negative curva- 
ture of the Bloch dispersion curve at the momentum mv/H. 
This is because in the reference frame where the lattice is at 
rest, particles' momenta are shifted by the value mv/h, and 
the effective mass, given by the second derivative of the curve 
in this point, becomes negative. This implies that the effective 
kinetic energy of particles and the interaction energy have dif- 
ferent signs, which immediately leads to the instability of the 
repulsive Bose gas at sufficiently large numbers of particles. 
For a shallow lattice, the Bloch dispersion curve is very close 
to that of free particles, except for momenta near the edge of 
the Brillouin zone k — zLkG/2, where the dispersion has a 
negative curvature (see, e.g., Ref. [72]). Hence, the condition 
for the dynamical instability takes the form fcc/S = mv/h. 
This coincides with the above condition for the break point, 
obtained within the generalized Landau criterion of superflu- 
idity, suggested in Sec. 111. A similar analysis can also be 
carried out for the TG gas. 

The frictionless motion at some values of the parameters 
V, 7, a is consistent with the presence of persistent currents 
in the ID Bose-Hubbard model [73-78]. As discussed in 
Sec. Ill, the drag force can be considered as a measure for su- 
perfluidity in the absence of the order paramenter Then Fig. 9 
represents the phase diagrams in the u-7 and v-l/a planes. 
They are similar to that of Ref. [75]. One can see from the 
diagrams that there is no sharp transition from superfluid to 
isolated phase in ID, which is consistent with the experimen- 
tal findings of Ref. [19]. Note that in paper [75], superfluidity 
was examined in terms of quantum phase slips, discussed in 
Sec. V E. So, the both quasiparticle and quantum phase slip 
descriptions lead to the same results. 



B. ID Bose gas in a moving random potential 

Recently, direct observation of Anderson localization with 
a controlled disorder was reported in a one-dimensional 
waveguide [20]. An initially trapped Bose gas is released into 
a one-dimensional waveguide, where it is exposed to an arti- 
ficially created speckled potential. The wave function of the 
spreading atoms can be directly measured by optical methods. 




FIG. 10. (Color online) The same diagram as shown in Fig. 9, but 
here the drag force is represented as a function of velocity and inverse 
filling factor. 



Anderson localization was initially predicted for a system 
of non-interacting particles [79] . Modern experimental tech- 
nique in cold atoms allows to control the strength of interpar- 
ticle interactions and the parameters of speckled potentials as 
well (see, e.g., reviews [80, 81]). Delocalization of a disor- 
dered bosonic system by repulsive interactions was observed 
[22]. Thus, the interplay between interactions and disorder 
remains a fundamental problem in this field. Considerable 
attention has recently been given to this problem in many pa- 
pers [82-91]. Most of the papers consider weakly-interacting 
bosons or the regime of infinite interactions (TG gas). By 
using the generalized Landau criterion of superfluidity, in- 
troduced in Sec. Ill), one can study the behaviour of a ID 
Bose gas in weak random potentials and obtain the superfluid- 
insulator phase diagram for arbitrary strength of boson inter- 
actions. 

It is convenient to consider random potentials by introduc- 
ing an ensemble of various potentials. Then the random poten- 
tial features can be obtained by averaging over the ensemble. 
One of the most important characteristics of random poten- 
tials, created with speckled laser beams, is their correlation 
function {Vi{x)Vi{x')) — g{x — x'), where (• • • ) stands for 
the average over the random potential ensemble. For an arbi- 
trary potential profile, the drag force is calculated with Eq. (8). 
Taking the average of this equation with respect to the ensem- 
ble, we obtain the drag force acting from the moving random 
potentials. At zero temperature, we derive 

/■2fcc 

(Fv) = / dkkg{k)S{k,kv). (39) 
Jo 
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FIG. 11. (Color online) Zero temperature phase diagram for 
superfluid-isolator transition of the Bose gas in a moving random po- 
tential: drag force (in units of 2-KmV^arN/h^) versus the potential 
velocity (in units vf) and the inverse correlation length (in units fcp). 
The dashed (blue) lines correspond to the intersection of the limiting 
dispersions ijj± and the transition line: v = aiiu±{2/ai)/2. On the 
left to the lines, the drag force is zero, because within the integration 
limits of Eq. (39), the values of the DSF is zero (see Fig. 4). 

Here g{k) = {\Vi{k)\'^) / L is the Fourier transform of the cor- 
relation function g{x). The integration in (39) is limited, be- 
cause for speckled laser beams, the function g{k) always has 
a finite support due to the limited aperture of the diffusion 
plate generating the random phase [92, 93]. So, g{k) — for 
\k\ > 2kc- For estimations, we take a realistic correlation 

function [92, 93] g{k) = TTVia,e(l - ^) (l - 

Here <Tr = 1/ kc is the random potential correlation length, 
depending of the parameters of the experimental device. It is 
usually of order of 1 /im [81]. Thus, we have three parame- 
ters governing the phase diagram: the potential velocity, the 
interaction strengh, and the coiTelation length. 

Like the lattice potential, the random potential is non-local, 
and, as a consequence, the drag force (39) is proportional to 
the total number of particles. The results obtained with the 
interpolating formula (26) are shown in Fig. 11. As for the 
shallow lattice potentials, the figure can be treated as the zero- 
temperature phase diagram for the superfluid-insulator tran- 
sion: superfluidity assumes zero or strongly suppressed values 
of the drag force. 

In the Bogoliubov and TG regimes, one can find the drag 
force analytically. For 7 ^ 1, we obtain from Eqs. (14) and 
(39) 

(Fv) = i^oe(c-ec)e(i-i)(i-<). (40) 



Here Fq = 2TrmV^(7rN /ft^ is a unit of force, t = 
-KTKj^yJi^ — Q., and = \fll''^ is the sound velocity (in 
units of Up) in the Bogoliubov regime. In the TG regime, 
Eqs. (11) and (39) yield 

(Fv> =Fo[/i + (/2-/i)e(A+-Ao)-/2e(A_-Ao)], (41) 

where we introduce the notations Ag = 2/(7rncri.), A± = 2|^± 

i|, A ^ - - h ^ 

The obtained results are in accordance with the existence 
of a mobility edge for a particle moving in a random poten- 
tial with a finite correlation length g^. In this case, the mo- 
bility edge is given by [83, 85] fc„iob = kc = l/ffr- If 
|fc| > kc, then the fc-wave propagation is not suppressed, 
while in the opposite case, the particle wave function is lo- 
calized (Anderson localization), which leads to the particle 
immobility. In the TG regime, the gas is equivalent to the 
ideal Fermi gas, where the mobile particles lie in the vicinity 
of the Fermi points. In the reference frame where the ran- 
dom potential is at rest, the absolute value of momentum of 
the lowest Fermi level is given by hk^ ~ \hk-p — mv\. When 
> ^c, the system should be superfluid. This is consis- 
tent with Eq. (39), which yields zero value of the drag force 
for /cp > kc- Indeed, the value 2A:p corresponds to the cross 
point of the dispersion curves for the TG gas and the transition 
line: w±(2A:p) = 2A:pti, where ijj± is given by Eq. (17) with 
infinite 7. Once 2fcp > 2kc then the value of drag force is 
zero, because the DSF is zero within the limits of integration 
in Eq. (39), see Fig. 4. 



VIII. CONCLUSION 

The absence of a well-defined order parameter makes the 
behavior of the ID Bose gas rather unusual in comparison to 
the 3D case. As shown in Sec. II, the ID Bose gas exhibits 
superfluid phenomena of equilibrium type (Hess-Fairbank ef- 
fect, analogous to the Meissner effect in superconductivity) 
but in general does not show dynamic superfluid phenomena, 
such as persistent currents in a ring. Instead of a phase tran- 
sition to full superfluidity as it is known in 3D, the ID Bose 
gas shows a smooth crossover and reaches the metastability 
of cuiTents only in a weakly-interacting limit. In this case, 
the drag force, being a simple integral parameter, can be cho- 
sen as a quantitative measure of superfluidity. Superfluidity 
assumes zero or strongly suppressed value of the drag force 
for the Bose gas moving in different small external potentials. 
Which value of the drag force should be taken as a the transi- 
tion threshold becomes a question of convention in ID. 

The drag force turns out to be fundamental, because it 
generalizes the Landau criterion of superfluidity. The gen- 
eralized Landau criterion, based only on energy and momen- 
tum conservation, does work when the usual Landau criterion 
fails. This is because the drag force effectively involves not 
only the spectrum but the probability of transitions to excited 
states. A good example is the dynamical instability of weakly- 
interacting ID bosons, moving in a shallow lattice. As shown 
in Sec. VII A, the generalized Landau criterion not only sue- 
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cessfully predicts the dynamical instability of the system but 
gives the quantitative characteristics of the phase transition. 

It should be noted that the suggested approach has an appar- 
ent disadvantage. Being based on first order time-dependent 
perturbation theory, the scheme cannot describe changes of 
the ground state caused by the perturbing potential. Despite 
of this fact, it can describe well the superfluid-insulator phase 
diagram, when the system propagates through shallow lattices 
or random potentials. 



ACKNOWLEDGMENTS 

The authors are grateful to Antony Leggett, Lev Pitaevskii, 
Hans-Peter Biichler, Audrey Kolovsky, Anatoly PoUcovnikov, 
and Andrew Sykes for fruitful discussions. 

J-SC acknowledges support from the FOM foundation of 
the Netherlands. AYC acknowledges a visiting grant from 
the ESfSTANS activity of the European Science Foundation 
and the support from the JINR-IFTN-HH projects and thanks 
Massey University for hospitality. JB was supported by the 
Marsden Fund Council (contract MAU0706) from Govern- 
ment funding, administered by the Royal Society of New 
Zealand. 



Performing the integration over y in the matrix element of 
the interaction Hamiltonian H' = Vi{\y — Xj\) yields for 

|(fef,n|H'|fei„,m)|' = \Vi{q)\^ \{n\Spg\m}f /V^ . (A4) 

At nonzero temperature, the medium can be in an arbitrary ini- 
tial state |m) with the statistical probability exp{—PEm)/2- 
In this case, the transition rate (Al) should be averaged over 
the statistical ensemble, and we arrive at 

'"^^^ = Y^v^^\''^-''-2^)- ^^^^ 

Here we use Eq. (A4) and the definition of the DSF (4). 

In order to obtain the energy loss per unit time, we need to 
sum up the energy transfer (A3) weighted with the rate (A5) 
over all the final states 

E = -Y^nj{q)(^hq-v-^). (A6) 

Replacing the sum by the integral in the thermodynamic 
limit J2q V/{2'7r)° J d^q and subsituting Eq. (A5) yield 
Eq. (3). 



Appendix A: General formula for the drag force from Fermi's 
golden rule 

Once an impurity moves into a homogeneous medium with 
velocity v, it is scattered by the medium particles. In general, 
the scattering leads to transitions with momentum and energy 
transfer and, consequently, to a finite value of energy loss per 
unit time. It can be calculated within linear response theory 
[5, 39, 94] or from Fermi's golden rule [94, 95]. 

During the scattering, the initial state of the composite sys- 
tem is assumed to be |fein,m) = |fein)|m), where |m) is 
initial state of the medium, the incident particle state is the 
plane wave |A;) = exp(ifc • y)/Vv with wave vector fcin, and 
V = is the D-dimensional volume. The wave vector is 
connected to the initial velocity by the relation v = hki^/mi. 
The same notations are adopted for the final state | fef , n) = 
I fef ) |n). The rate for the scattering process (that is, the transi- 
tion probability per unit time) is given in the lowest order in 
the impurity interaction by Fermi's golden rule 



(fcin, fef) =^ ^ |(fef,n|i?'|fei„,m)|^ 

n 

X 5{En + Tk, -Em- TkJ 



(Al) 



where = h'^k^/{2m) is the energy dispersion for the im- 
purity. The momentum and energy transfer are given by 

hq = hkin — hkf, (A2) 
hu: = Tk,^-Tk,=nq-v-h^q^/{2mi), (A3) 

respectively. 



Appendix B: The drag force in the RPA approximation 
1. The general expression 

In this appendix, we will use the dimensionless variables 

A = q/k-p, V = hw/e-p, C = f^/'^F, where kp = im, ep = 
h^kp/{2ni), and vp = hk-p/m are the Fermi wave vector, 
energy, and velocity, respectively. Besides, it is convenient to 
introduce the small parameter a = 2/7, which can take the 
values ^ a ^ 1/4 within the RPA scheme [42]. In terms of 
the new variables, Eq. (9) reads 



/v 



+ 00 



dAAs(A,2^A), 



(Bl) 



where we put s(A, v) = £pS{kp\, £pv/h)/N . 

The DSF in the RPA approximation was calculated and de- 
scribed in details in Sec. IVB of paper [42]. The result at zero 
temperature can be written in the form 

S(A, U) = Sreg(A, u) + A{\)5{v - vo{X)) (B2) 
with the regular part of the dimensionless DSF 

A (l-3a)2(l-2a) 



4 [(1 - 3a)2A -ahlnf] + [a7r/i]2 



(B3) 



This part is localized in the same region | A^ — 2A1 (1 — 2a) ^ 
^ (A^ + 2A)(1 — 2a), as the DSF in the linear approxima- 
tion, given by Eqs. (16) and (17). Here we put by definition 



i/2-(A2-2A)2(l-2a)2 



(A2 + 2A)2(l-2a)2 



(B4) 
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and 



a 3a(l-8a/3) 



MA,.).(1-W4)----^^, (l-2«)^ 



(B5) 



The quantities in Eq. (B2) are given by A{X) = A{k-p\)/N, 
1^0 (A) = lThjQ{k-pX) / e-p with A{q) and wo('z) being defined 
exactly in Sec. IVB of Ref. [42]. By employing the general 
formula (Bl) and the RPA expression (B2), we can evaluate 
the drag force by numerical integration. The 5-function con- 
tribution appears as a result of intersection of the curve t'o(A) 
and the line 2^A in the v-\ plane and does not play a role at 
small values of the impurity velocity, see Fig. 12. One can 
easily show that the (5-contribution appears for a ^ 2/9 when 
s; C < ore ^ where 



r(a) = (l-2a). 



(1 - 3a)Va2 - 1 



3(l-8a/3) + (l-3a)2/a'' 



C±(a) = (l-2a) 



4 - 9a ± V^Vie - 71a + SOa^ 
4(l-3a + 2a2) ' 




0.00 0.05 0.10 0.15 0.20 0.25 



FIG. 12. (Color online) The dark (blue) region shows where the 5- 
function contribution exists. Such a contribution results from the 
intersection of the curve uo{X) and the line 2^A in the v-\ plane. 
One can see that there is no (5-function contribution at sufficiently 
small velocity of the impurity. Here = c/uf, and — 1 — 2a in 
the strong-coupling regime. 



2. The drag force for small velocities at zero temperature 

Since at the small velocities the linear approximation (18) 
fails, we need to apply here the full RPA expression for the 
drag force. At small values of ^ = v/v-p, the limits of the inte- 
gral in Eq. (Bl) becomes very close to each other, and we can 

put f = r2(i+«/«=) 

and obtain from Eq. (B3) 



'AXXs{\v = 2eA) ~ srcg(2, 408^/^0 



4^(1 -3a)2 



Z32 



5a 



e 3a 



2 V-^ 3 / 



2 ' 



(B6) 



where 



D = 2(l-3a)^ 



2-5a- 



In 



(4-')- 



^^30/ 8aN 

In principle, the logarithmic term should dominate in the 
denominator when ^ ^ and one can neglect all the other 
terms; however, this approximation works only at very small 
values of S^. A much better approximation can be obtained 
by keeping the logarithmic term together with the zero-order 
terms in ^: 



/v 



(1 - 3a)2 



(1 - 5a/2)2 



(l-3a)^ 
(1-5q/2) 



2aln 



(B7) 
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